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Abstract 

oo : 

■ Starting with the vacuum fluctuation, it is known that gravitinos will be created 

just after inflation, with number density ~ 10 -2 M 3 where M is the mass of the 
inflaton. Here, we argue that creation may be expected to continue, maintaining 
about the same number density, until a usually much later epoch. This epoch is 
either the 'intermediate epoch' when Hubble parameter falls below the gravitino 
mass, or the reheat epoch if that is earlier. We verify that such late-time creation 
indeed occurs if only a single chiral super field is relevant, using the description of 
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. the helicity 1/2 gravitino provided recently by Kallosh et. al. ( hep-th/9907124| ) 

^5 



and Giudice et. al. ( liep-ph/9907510|) . Arguments are presented in favor of late- 



time creation in the general case. For the usual inflation models, M is rather large 



and gravitinos from late-time creation are so abundant that a subsequent era of 
(U '. thermal inflation is needed to dilute them. 



Introduction Gravitinos are created in the early Universe by thermal collisions 0, and 
some time ago |2| it was pointed out that they will also be created non-thermally, starting 
from the vacuum fluctuation that exists well before horizon exit during inflation. To see 
whether creation from the vacuum is significant, one needs equations which describe the 
evolution of the gravitino mode functions. Such equations have recently been presented 
II, |j], for the case that only one chiral superfield is relevant. Using them, the authors find 
that the number density of gravitinos created just after inflation is of order 10 _2 M 3 , where 
M is the mass of the inflaton after inflation, and they have noted that these gravitinos may 
be more abundant after reheating than those created from thermal collisions. A similar 
result has been obtained for hybrid inflation, which involves at least two fields, assuming 
that in this context the gravitino becomes the goldstino of global supersymmetry || (see 
also ||). 

In this note, we argue that this is unlikely to be the end of the story J7[. Rather, 
creation is likely to continue, maintaining about the same number density, until either 
the 'intermediate epoch' when Hubble parameter falls below the gravitino mass, or the 
reheat epoch if that occurs earlier. We begin by verifying that such late-time creation 
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indeed occurs if only a single chiral superfield is relevant, using the the description of 
the helicity 1/2 gravitino provided recently in P, |[]. Then we consider the case where 
other fields are relevant; a different field to break supersymmetry in the vacuum, a third 
field to allow hybrid inflation, and the fields corresponding to particles created just after 
inflation by preheating. In all these cases, we argue that late-time creation will continue 
until the intermediate epoch, unless it is terminated by preheating. Finally, we consider 
the cosmological consequences of late-time creation, which are rather dramatic in the 
most popular models of inflation. 

Describing the gravitino in the early Universe To calculate the abundance of 
gravitinos created from the vacuum, one needs equations describing the evolution of the 
mode functions for the helicity 1/2 and 3/2 gravitino states, as seen by a comoving 
observer in the expanding Universe. 

For the helicity 3/2 state, one can safely proceed by using the Rarita-Schwinger 
equation with appropriate constraints, evaluated in Robertson- Walker spacetime, and 
with time-dependent mass m 3 / 2 (t) given by the usual N — 1 supergravity formula. (We 
denote the vacuum value of m 3 / 2 (£) by m 3 / 2 without an argument.) A suitably-defined 
mode function satisfies the spin 1/2 equation, with the same effective mass m 3 / 2 (t) that 
appears in the Rarita-Schwinger equation || [|, |], [|. As m 3 / 2 (t) can hardly be bigger than 
the Hubble parameter H if it is to vary non-adiabatically,[| the abundance of helicity 3/2 
gravitinos created from the vacuum will no bigger than that of modulini created from the 
same mechanism, which is negligible [§, |n| compared with the abundance of gravitinos 
from thermal collisions. 

For the helicity 1/2 state, the same procedure may in general require modification, 
because the would-be goldstino (existing in the limit of global supersymmetry) may be 
a time-dependent mixture of the spin 1/2 fields. To avoid this problem, the case where 
only a single chiral supermultiplet is relevant has been studied |3|, [|. It is found that a 
suitably-defined mode function again satisfies the spin 1/2 equation, but with a different 
effective mass rh(t). 

In a general model of inflation, the idealization of a single chiral supermultiplet will 
not be adequate. In hybrid inflation models one needs, in addition to the slowly rolling 
inflaton field, one or more additional fields to provide the constant part of the potential. 
Even if the inflaton field is the only one relevant for inflation, a second field is generally 
needed to break supersymmetry in the vacuum. Even so, the possibility that the inflaton 
field is the only relevant one is not excluded. Also, the field responsible for supersymmetry 
breaking may at some stage oscillate and dominate the energy density of the Universe, 
even if it is not responsible for inflation. Let us proceed on the assumption that only a 
single field is involved in gravitino creation, so that the equations presented in [EL [| can 
be used. 

lr The N — 1 supergravity expression for H 2 contains a term — m^ 2 (t), which can hardly vary rapidly 
if it is canceled to high accuracy by the rest of the expression. Of course, a precise cancellation does 
occur at the present time, for some unknown reason (the cosmological constant problem). 
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Gravitino creation from the oscillation of the supersymmetry-breaking field 

The effective mass m(t) of the helicity 1/2 mode function is given by P, |J 2 

fh(t) = m V2{t) - ^m 3 / 2 (t)(l + M) - hiA 2 -n, (1) 

where 

= P-3M 2 m 2 /2 (t) 

1 " p + 3M 2 m 2 /2 (t) [ } 

= 2 3M P m 3 / 2 (t) 

2 " 3 p + 3M 2 m 2 /2 (t) K ' 

A = A 1 + iA 2 = e ix (4) 
1 

» = 2 X ' w 

In this expression, Mp = 2.4 x 10 18 GeV is the Planck scale. 

The energy density p and the pressure P are supposed to be dominated by a real 
scalar field <fi. Taking its kinetic term to be canonical, 

p = V + ±<j> 2 (6) 

P = -V+\P, (7) 

where V{4>) is the potential. The continuity equation is p = —3H(p + P), equivalent to 
the field equation 

<j> + 3H<j) + V = (8) 

where H is the Hubble parameter, related to the energy density by 3MpH 2 = p. 

For future reference, let us note that the equations for A±, A 2 and p may be written 
in terms of H, w = P/p, and m 3 / 2 (t), with the time- derivative of this last quantity 
eliminated using Eq. (H). One finds 



wH 2 - m 2 3/2 (t) 
H 2 + ml /2 (t) 



2 In Eq. (3.23) of |3j, the signs of the second and third terms should be reversed jll), after which that 
equation becomes identical with Eq. (|l|) upon making the identification a^ 1 ^!^ = —m(t). The sign of the 
effective mass is not physically significant, and in particular it will not affect the gravitino abundance. 



(Equivalently, the sign of the last term of Eq. (12) is not physically significant, since it corresponds to 
replacing a mode function by its complex conjugate.) We have checked that Eq. (|l|) is equivalent 
to the more complicated expressions given in [Q. In both || and Q], it is assumed that the field is 
canonically normalized, but the results given there are valid also for arbitrary normalization [^2[ . In 
fact, for a single real scalar field, one can always transform at least locally from arbitrary normalization 
to canonical normalization. In the supergravity model, (f> is the real part of a complex field, and the 
corresponding transformation of the complex field is at least locally holomorphic, leading to an equivalent 
supergravity theory. 
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w 2 + 2(1 +w) m 2 /2 (t)/H 2 



l + m 2 (t)/H 2 



(10) 



"3/2 V 

3 /W1 . N lw -3F(1 + w)(w -At) 

For momentum k/a, a suitably defined helicity 1/2 mode function satisfies the equa- 
tion 

u" + (k 2 + (am) 2 + i(am)'^) u = , (12) 

where the prime denotes differentiation with respect to conformal time drj = dt/a, and 
a is the scale factor of the Universe such that H = a/a. To calculate the gravitino 
abundance, one starts at early times with the negative-frequency solution u oc e~ luJV , 
corresponding to the vacuum. At late times there is a linear combination of positive and 
negative frequency modes, and the occupation number is the coefficient of the positive 
frequency mode. Significant production occurs with momentum k/a if there is appreciable 



violation of a weak adiabaticity condition JTD 



| (am)' | < or = Ar + (am) , (13) 

where the average is over a conformal time interval uj~ 1 . In practice fc max , the biggest k 
for which significant creation occurs, is simply the biggest value achieved by am, within 
the regime where m varies non-adiabatically (|m| > m 2 ). 

Let us follow the evolution of m, to estimate fc max - During inflation, m ~ m^^it) 
is much less than H in magnitude. After inflation, oscillates about its vacuum value 
which we set equal to zero so that <fi ~ (f>o(t) sinMt where M 3> H is the inflaton mass. 
We now have m ~ — /x = — with 

|0 2 -V-3Mgm 2 /2 (t) 

COS Y = — rr (14) 

\cj) 2 + V + 3M 2 m 2 3/2 {t) V ; 

cos 2Mt - m 2 ,Jt)/H 2 

l + m 2 /2 (t)/H 2 ■ 1 ] 

When = 0, cosx = — 1- The crucial point is that the maximum value of cosx is less 
than +1, because cosx = 1 would correspond to V — — 3Mpm 2 / , 2 (t) < 0. It follows that x 
oscillates in some range e < x < 27r — e. With the reasonable assumption m 2 / 2 (£) ~ m 2 / 2 , 
e ~ 1713/2/ H 1. While x is rising, m ~ —M, and while it is falling, m ~ +M. Because 
of the sign switches, the evolution of the mode function is non-adiabatic, and gravitino 
creation continues, with every- increasing comoving momentum k ~ aM. 

When H becomes of order m 3 / 2 , the sudden switches in sign of m give way to a 
complicated variation, still on the timescale M _1 and with typical value of order M. (If 
desired, this variation is conveniently calculated from Eqs. ©-(0), with w = cos2Mt.) 
Creation finally ceases only when H falls below m 3 / 2 , and according to Eqs. ([I]) and ( |TT] ) 
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in falls smoothly to the true value 7713/2, recovering the flat spacetime description of the 
gravitino. 

We studied this example primarily to provide an existence proof, that late-time grav- 
itino creation occurs in at least one case. This is the case that the oscillating field 
responsible for the energy density of the Universe, is the same as the field breaking su- 
persymmetry in the vacuum. Before moving on, we note that this case may by realized 
in Nature if supersymmetry breaking is gravity mediated, though the mass M would 
then probably be too small for late-time creation to be significant. Indeed, superstring 
theory suggests the existence of several fields (moduli) with only gravitational strength 
interactions, and mass very roughly of order m 3/ / 2 . (Kahler stabilization of the moduli 
can give masses several orders of magnitude bigger, for instance a dilaton mass of order 
10 6 GeV is found in JTBJ.) At least some of the moduli may well be displaced from the 
vacuum at the end of inflation, and one of them might be the inflaton). Any displaced 
modulus at first tracks the changing vacuum value, but it starts to oscillate at the inter- 
mediate epoch [n| , and then dominates the energy density. Finally, one of the moduli is 
usually supposed to be responsible for supersymmetry breaking in the gravity-mediated 
case. The case that we have studied will be realized if this modulus is also the inflaton, 
and is the only one that oscillates. As far as the oscillation regime is concerned, the case 
is also realized if the inflaton is different from the supersymmetry breaking modulus, but 
decays early so that the latter is the only oscillating field at the intermediate epoch. 



Late-time gravitino creation in the general case When additional fields and/or 
particles are involved, equations describing the evolution of the gravitino field are not 
yet available. However, late-time creation seems likely to occur quite generally up to the 
intermediate era, unless reheating occurs first. 

Consider first models where the slowly-rolling inflaton field is the only one relevant 
for inflation, but something else breaks supersymmetry in the vacuum while ensuring 
that the potential vanishes there. In the model where the inflaton field did both jobs, 
late-time creation occurred because supergravity corrections to global supersymmetry 
become very important every time the potential the potential dips to zero, causing the 
mode function to vary non-adiabatically. One should not expect that the introduction 
of something else to break supersymmetry, would restore the adiabaticity. 

Next, consider hybrid inflation models, where two fields are oscillating after inflation, 
with presumably a third breaking supersymmetry in the vacuum. In supersymmetric 
models, the two oscillating fields typically have the same mass M after inflation. In gen- 
eral they are not oscillating in phase (though see Jl5| for a case where they are), which 
means that they cannot be replaced by a single field. Just after inflation, gravitino cre- 
ation in this situation may be estimated in the global supersymmetry limit, as described 
in || for the case of F-term models. The number density will again be of order M 3 , since 
the inflaton mass climbs to that value in a time of order M~ l . Again, there is no reason 
to expect gravitino production to stop, since supergravity corrections will again become 
important every time the potential dips to zero. 
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In both of these situations, a definite calculation will clearly become possible when 
the supergravity formalism of |3], [| is extended to include two or more chiral superfields. 
It is not so clear how to calculate things when particles as opposed to homogeneous fields 
become important, but still it is fairly clear what will happen. Consider first the case that 



preheating [16] converts most of the oscillating field energy into marginally relativistic 
particles. The energy of the latter will be somewhat reduced by redshifting even if it does 
not decay promptly into radiation, and one expects that after a few Hubble times oscil- 
lating fields again account for a non-negligible proportion of the energy density. Then 
late-time gravitino creation will presumably continue until the intermediate era, unless 
reheating happens first. After reheating, practically all of the energy density is in radia- 
tion, and one expects creation to stop, because nothing is varying rapidly. We emphasize 
that, although these expectations look quite reasonable, it is at the moment totally un- 
clear how to describe the gravitino, when supersymmetry breaking comes mainly from 
the particle gas in the early Universe. 

Cosmological significance of late-time gravitino production After production 
stops, the occupation number will be of order 1 below k = fc max , giving number density 



n ~ — -cT 3 / k 2 dk (16) 

47T 2 Jo 

~ 10- 2 (W«) 3 - (17) 

The number density when creation stops is n ~ 10 _2 p 3 , where p = fc max / a is the maximum 
momentum of the created gravitinos. 

We focus on the extreme case, where gravitino creation ends only at the in termediate 
era. This corresponds to energy density of order M|, where Ms ~ ^jMpm 3 / 2 is the scale 
of supersymmetry breaking. In gravity-mediated models of supersymmetry breaking, 
7713/2 ~ 100 GeV and Ms ~ 10 10 GeV. In typical gauge-mediated models, m 3 / 2 ~ 100 keV 
and M s ~ 10 7 GeV. 

At the intermediate era, the gravitino number density is of order 10~ 2 M 3 . The relative 
abundance at nucleosynthesis is therefore 



10-^. (19) 



Here, s is the entropy density at nucleosynthesis, and 7 _1 is the increase in entropy per 
comoving volume (if any), between reheating at temperature Tr and nucleosynthesis. If 
the entropy increase comes only from a late-decaying particle, there can be only a modest 
increase corresponding to 7 ~ Tfr/Teq, where EQ is the epoch when the particle first 
dominates the energy density, and FR is the final reheat epoch when it decays. This gives 
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> Tfr > 10 MeV, where the bound comes from nucleosynthesis. However, N e-folds 



of thermal inflation p7| , [bj, [L8[ may also occur. This would reduce 7 by an additional 
factor e~ 3N . One bout of thermal inflation typically gives M, |18| iV ~ 10 and a total 7 



perhaps of order 10 -15 . Two bouts might give N ~ 20 and up to 7 ~ 10~ 30 |L4 . 

The cosmo logical significance ]I[ of the gravitino depends on its true mass m 3 / 2 . With 
gravity-mediated supersymmetry breaking, one expects m 3 / 2 ~ 100 GeV to 1 TeV, and 
observation then requires 

n/s<10~ 13 . (20) 
The abundance of gravitinos from thermal collisions is then 

n/s ~ 10- 13 ( 7 T R /10 9 GeV) , (21) 

leading to the bound 7TR < 10 9 GeV.f] Using instead Eq. (|19"D, we find 

10-!! „ f^L_) 3 flOO^V < l. (22) 

s \WGeVj U0 9 GeVy^ m 3/2 J ~ y ' 

Gravitinos created from the vacuum are more abundant than those from thermal colli- 
sions, if M > 10 7 GeV. Without thermal inflation, we need 7TR > 10 MeV and therefore 
M < 10 11 GeV. In the worst case Tr ~ Ms, even a single bout of thermal inflation allows 
only M < 10 11 GeV. With two bouts, or with one bout and low Tr, more or less any M 
up to Mp might be accommodated. 

With gauge-mediated supersymmetry breaking, one expects very roughly 1 keV < 
m 3/2 ^ 100 GeV, with the upper decades disfavoured. Unless m 3 / 2 > 100 MeV, this leads 
to a stable gravitino with present density 

fi 3/2 - 10 5 (^Tr) ~ & 1 ( 23 ) 

3/ VIOOkeV/ s K J 

If creation from thermal collisions dominates, then very roughly [IJ 

3/2 I m 3/2 ) Uo 4 GeV7 1 } 



Using instead Eq. ( TT9| ) 



In this case, creation from the vacuum is more efficient if M > 10 8 GeV, and we need 
thermal inflation if M < 10 10 GeV. In the favored case m 3 / 2 ~ 100 keV, corresponding to 
Ms ~ 10 7 GeV, the constraints on M are similar to those in the gravity-mediated case. 



We refer here to the thermal creation at the initial reheating. If thermal inflation subsequently 



occurs, the most significant thermal creation occurs at the final reheating, and 72r in Eq. (21) is to be 
replaced by the final reheat temperature. 
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So much for the case that reheating happens after the intermediate epoch correspond- 
ing to Tr < Ms. By way of contrast, consider the opposite extreme of instant reheating, 
T R ~ V l/A . Then, 

n - ~ 10-7 (£)' • (26) 
s W*/ 

Either M/V 1 ^ should be sufficiently small, or thermal inflation is again needed. 

Conclusion Let us end by considering the implications of late-time gravitino creation 
for models of the early Universe, and in particular of inflation. Supersymmetric models 
of inflation are reviewed for example in fli~9"| |. In some of them, notably those using only 



flat directions, the mass of the field(s) oscillating after inflation is small, and late-time 
creation cannot be significant. However, in the most models, normalized to give the 
correct prediction for large scale structure, the inflaton mass is bigger than 10 10 GeV. 
Using the popular names, examples include chaotic inflation, the most popular hybrid 
inflation models (D term, and the usual F term model) and most new inflation potentials. 
The popular hybrid inflation models require M ~ V 1 ^ ~ 10 15 GeV or so. 

In all of these cases, the possibility of late-time gravitino creation drastically changes 
ones view about the reheat temperature. Starting the discussion with instant reheat, 



Eq. ( p6]) represents a significant constraint. The constraint becomes stronger as the reheat 
temperature is lowered, because gravitino creation persists to a later epoch so that n/s 
increases. Only when the reheat temperature falls below Ms, does the constraint, now 



represented by Eqs. fl22|) and (|25|), start to become weaker. In many cases, the constraint 
cannot be met simply by lowering the reheat temperature. Instead, the entropy dilution 
factor 7 has to be small, often so small that thermal inflation is required. 

To avoid these constraints, one might turn to models using only flat directions, leading 
to M perhaps of order and giving negligible gravitino creation. Such models include 
the modular inflation mentioned earlier (so far without a concrete realization in the 
context of string theory) and certain hybrid inflation models. 

The next step will be to check that late-time gravitino production actually occurs 
in the models mentioned, using equations that describe the gravitino in the presence of 
several fields. After that would come the more challenging task of describing the gravitino 
when a particle gas is the dominant source of supersymmetry breaking. 
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